Two electron entanglement in quasi-one dimensional system: Role of resonances 
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We analyze the role of resonances in two-fermion entanglement production for a quasi one- 
dimensional two channel scattering problem. We solve exactly for the problem of a two-fermion 
antisymmetric product state scattering off a double delta well potential. It is shown that the two- 
particle concurrence of the post-selected state has an oscillatory behavior where the concurrence 
vanishes at the values of momenta for virtual bound states in the double well. These concurrence 
zeros are interpreted in terms of the uncertainty in the knowledge of the state of the one par- 
ticle subspace reduced one particle density matrix. Our results suggest manipulation of fermion 
entanglement production through the resonance structure of quantum dots. 

PACS numbers: 03.67.Mn,03.65.Ud, 73.50.-h 



Entanglement production and quantification has been 
given much recent attention due to its importance 
as a resource for quantum information and quantum 
communication [l], In this direction, there have been 
recent proposals for producing bipartite fermionic entan- 
gled states in the solid state environment focusing on 
the role of direct interaction between particles. Some of 
these approaches involve direct coulomb interactions in 
quantum dots [3] and interference effects H, phonon me- 
diated interactions in superconductors fs", and Kondo- 
like scattering of conduction electrons [7]. Nevertheless, 
it has been shown that fermion entanglement can be 
achieved in absence of such interactions [8] in the form 
of particle hole entanglement even when fermions are in- 
jected from thermal reservoirs. In such a setup the orbital 
degree of freedom is entangled. Other implementations 
based on the non-interacting scheme have been proposed 
that entangle the spin degree of freedom and are thus 
more robust to decoherence [9] because of the weaker 
coupling of the spin to the environment. 

In this work we address the problem of entanglement 
generation for electrons in the context of a two channel 
quasi-one dimensional conductor [l^ , following the scat- 
tering matrix formalism of reference [8] . For the scatter- 
ing region, we choose a double delta potential, separated 
a distance d. Such a potential is the simplest, that ex- 
hibits resonances, that can be analyticall y h andled. The 
problem is solved for the concurrence [Ul m exactly, for 
all values of the barrier heights and separation as a func- 
tion of the incoming electron momenta. The concurrence 
of the entangled post selected state is found to oscillate 
while its envelope decays as a function of electron mo- 
mentum {ki) difference A/c = /c2 — ki. We find that the 
concurrence is exactly zero when one or both of the k val- 
ues hits the resonant states for the potential well. The 
concurrence zeros are then interpreted in terms of the 
uncertainty of the state in the one particle subspace by 
obtaining the reduced density matrix. We thus determine 
the role of resonances on the entangling properties of the 
well demonstrating new possibilities for fermion entan- 
glement control. We consider in detail the independent 



channel scenario but quantitative changes due to chan- 
nel mixing will be briefly discussed. Here we ignore the 
effects of temperature since they have been assessed in a 
general way in reference [13] and will not change qualita- 
tively the results reported here if a critical temperature 
is not reached. The set up for the system considered is 
depicted in Fig. [1] where a two electron wavefunction is 
injected at the left of a quasi-one dimensional conductor 
wire. The electrons move freely until they enter the in- 
teracting domain with potential V{x^y). Each electron 
is considered to pertain to a separate channel in the in- 
coming lead and gets transmitted or reflected within the 
same two channels. 

Let us set up the problem in a first quantized descrip- 
tion. Schrodinger's equation is given by 



2m \dx'^ dy'^ 



V{x,y) 



tlj{x,y) = Eilj{x,y). 



(1) 

The potential V{x, y) acts in a finite region of the coor- 
dinate X (see Fig. [1]). The boundary conditions on the 
wire are such that ?/^(x,0) = tlj{x^w) = 0. In the free 
regions (the leads) V{x^ y) = and using the definitions 
k'^ = 2mE/h^ and U (x, y) = 2mV(x, y)/fi^^ we can write 
the Schrodinger equation as 



dx'^ dy'^ 



^{E\\,x)xn{y)=^. (2) 



where /c^ = /cy + ^1 n- eigenfunction (/)(£^||,x) is 

given by e^^^w^ / ^2T:h?k\\/m^ ^±,n = rm/w and Xn{y) = 
^J^^inK^^ny- The integer n denotes the channel num- 
ber. When /c^ > ^ both channels are open. 

In the potential region i){x,y) = En=i V^n(^)Xn(y) 
where, if we define Umni^^) = Jo Xn{y)U{x,y)xm{y)dy, 
we get the system of coupled equations 



dx'^ 



V^n(^) = ^ Umn{x)^lJm{x). (3) 
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We denote the difference k'^ — K\ ^ = /c^ , omitting the 
suffix II, and we always understand that the difference is 
positive. 

We now fix n = 1,2 and choose Umn{x) = Umn'^{x)^ 
with v{x) = 5{x — d/2) + d{x + d/2) and Umn = "Kim- 
To begin with, we take Ui2 — U21 = which means no 
channel mixing. The composition of two delta scatterers 
at X = —d/2 and x = d/2^ in series, with corresponding 
scattering matrices Sj and Su, gives the symmetric 5* 
matrix 



(4) 



where 1 is the 2x2 identity matrix. 

We can now use Beenakker's[8j approach to arrive at 
the expression for the output wavefunction. Using the 
same notation 




aL,i(^)"ln,2(e)|0), 



(5) 



where aj^ ^ creates one electron on the left in channel 

/. Now bl^j{e) creates an electron in channel j incident 
from the right so that in matrix notation the input state 
can be written as 



|*m) = 



b. 



I^y 




|0), 



where the vectors are 4x1 and the matrix is 4 x 4, because 
there are two channel indices on the right and the left. 
The relation between input and output channels is given 
by the scattering matrix 



^out 
bout 



r t' 
t r' 



The entries r, t, r' and t' are 2x2 refiection and trans- 
mission matrices. After some algebra one arrives at the 
exact relation 

K^^]l2ClC2)|0). (6) 

For no channel mixing, and in terms of our particular 
potential, the r and t matrices are given through 



2^2iknd 



-iknd 



^3 



' 33 



1 



""33 



1 



and ri2 = r2i = ti2 = ^21 = 0, where the indices of 
the reflection and transmission amplitudes refer to the 
channel, and r^, tj (j = 1,2) have the expressions Vj = 
(ujj/2ikj)/(l-Ujj/2ikj) and tj = 1/ (I - Ujj /2ikj). 
The new element here is that now we have energy de- 
pendent transmission and reflection amplitudes and the 
presence of resonances because of virtual states in the 
barrier. 



¥.=<i<2|0) 



w 




V(x,y) 

FIG. 1: The scattering setup for a quasi-one dimensional wire 
of width w and a scattering region with potential V{x^y) con- 
sisting of a sequence of two delta potentials separated by a 
distance d . A two electron antisymmetrized wave function is 
injected at the left. The outgoing products according to Eq|6] 
consist of three terms: a) two electrons are reflected b) two 
electrons are transmitted and c) one electron is transmitted 
and the other reflected. 



In order to derive from the scattering result of EqlBlan 
entangled state one must now post-select or project out 
the appropriate component. In this case one can post- 
select by coincidence measurements where electrons are 
detected simultaneously at opposite branches of the dou- 
ble barrier, a well known experimental tool in optics [l3|. 
The useful term is flrst order in t and r generating par- 
ticles on both sides of the double barrier 



i$) = 



aL7&L|0>, 



(7) 



where 7 = rcjyt^ . The state is appropriately nor- 
malized. In order to compute the concurrence we use 
the convenient deflnition[15], that reduces the problem 
to identifying the W matrix in the expansion \^) = 
Za,f3'^o.p aibl\0), where a,p e {1,2,3,4} and Wap 
can be assumed antisymmetric. Expanding the product 
«Lt7^L|0) one flnds 



721<^Lt,2^Lt,l 



7l2<^L,lCt,2 + 722aL,2C^,2 

Then the antisymmetric part of Wa(3 is given by 



W 



2VTr77t 



/ 711 712 

721 722 

-711 -721 

V -712 -722 



The expression for the concurrence is rj = |(^^|^)| = 
£:"^^^Wa/3W^jy, where e^f^^^ is the totally antisymmetric 
unit tensor in 4 dimensions. Then 77 = 8IW12W34 + 
Wi3 W42 + Wi4 W23 1 . Computing r] for the W^ matrix 
above gives 2|det7|/Tr77'''. For the general case including 
channel mixing the matrix 7 is given by 



^12^11 
^22^11 



■ ^11^12 ^12^21 
^21^12 ^22^21 



^11^22 
^21^22 



(8) 
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FIG. 2: Concurrence as a function of the wavenumber differ- 
ence Ak = k2 — ki (in 27v/d units) without channel mixing. 
The barrier heights have been fixed to uq = {27v/d) • (1/100) 
and /ci,/c2 take values as shown. In the bottom panel the 
concurrence is zero whenever k2 hits a resonance while ki is 
in between resonances as indicated in the inset. In the top 
panel now ki is in the vicinity of a resonance and a concur- 
rence maximum occurs when a k2 closes onto another (going 
exactly to zero when the resonance is hit). 
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FIG. 3: (Color online only) Concurrence as a function of both 
ki and Ak — k2 — ki. The strength of the delta potential has 
been fixed dX uq — 2n /d (1/100). Limiting behaviors depicted 
in figure 2 connect smoothly as ki is varied. The zeros of 
concurrence always correspond to k2 hitting a resonance. 



It is obvious that when there is no channel mixing the 7 
matrix is anti-diagonal and the resulting concurrence is 
then 



2|r22||tii||rii||t22| 
|r22P|tiiP + |rnP|t22P 



(9) 



Substituting the values for the reflection and transmis- 
sion amplitudes one can derive 77 as a function of the 
incoming wavevectors. Since both transmission and re- 



flection matrices are /c-dependent, the resonant proper- 
ties of the device are relevant for the two particle concur- 
rence. It is worth noting a crucial point: If one computes 
de concurrence without post-selecting the result is zero. 
This is consistent with the fact that local processes can- 
not change the entanglement, which for the input state 
is zero. 

In the absence of channel mixing we depict, in Fig. 
[21 7^ as a function of the difference in the magnitude of 
wavevectors Ak = k2 — ki. We have scaled all wavevec- 
tors and the magnitude of the delta potential so they are 
in units of 27r/d. Having fixed the height of the barriers 
and the distance between them, the one particle reso- 
nances occur at fixed values shown in the figure panels 
as insets in units of 27T/d. Two well defined limiting be- 
haviors occur: The bottom panel depicts the case where 
ki = 27r/d; here we notice that as a function of A/c, r] 
shows an oscillating and decreasing pattern. The min- 
ima, which are exactly zeros of the concurrence, occur 
when k2 values hit a resonance while ki is in between 
resonances as indicated by the inset. Succesive zeros co- 
incide with succesive one particle resonances. 

The top panel of Fig|2l shows a second behavior occur- 
ing when ki is in the vicinity of a resonance, as indicated 
by the panel inset. The concurrence then is only appre- 
ciable within the resonance width, and goes to zero, ex- 
actly, when /c2 wavevector hits resonances. The full range 
of behaviors described and their crossovers between that 
of the bottom panel and top panel in Fig|2l are shown in 
Fig|3]in a representative range of ki and Ak. 

A useful tool to gain intuition on entanglement is to 
obtain the reduced one particle density matrix of the 
state in Eq|71 If there is entanglement, the resulting 
density matrix represents a mixed state showing there 
is uncertainty in the state of the particle. Vanishing of 
entanglement is then evidenced by the certainty of a par- 
ticular state. We stress though, that there is no new 
information regarding entanglement that is not already 
assessed in EqlQl Setting up the two particle density ma- 
trix p2 = |^)(^| and the tracing over one of the particles 
results in the matrix 



pi Trip2 



p2 rp2 

^22-'- 11 

2Tr77t 

TD2 rp2 

2Tr77t 



|7^l)(7^l| + |/:l)(/:l|) + 
|7^2)(7^2| + |/:2)(/:2|) (lo) 



where Rf^ = Ir^ip and analogously for Tu. The kets are 
defined as bl^^ -lO) = |7^i) and a^^^ .|0) = This is a 
mixed state (as can be seen by tracing over pf) where 
the remaining electron is projected, with probability 
^22^11 /2Tr77t, onto channel 1 on the right, |7^l). Note 
that for arbitrary reflection and transmission probabili- 
ties, the electron can be in any of the two channel states 
(1 or 2) signaling entanglement between the two elec- 
trons. This indeterminacy is destroyed once we hit a sin- 
gle particle resonance so that Tu = 1 (so that Ru = 0) 
becoming a certainty since pi = 1 /2 ( | T^i ) (JZi \-\-\jCi){Ci\) 
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and the electron can only be in channel 1. This explains 
the zeros of the concurrence at the single particle reso- 
nances. 

Some additional features of the figure can be accounted 
for using the above expression: The first maximum 
{ki = k2) in the bottom panel of Fig|2] corresponds to 
the maximum uncertainty in distinguishing one particle 
from the other. In this situation, the probability ampli- 
tudes for transmission (refiection) t (r) through either 
channel are the same (see Eq. [9]). As the wavevector dif- 
ference increases the concurrence envelope function drops 
monotonously indicating the uncertainty is also reduced. 
This can be seen from EqUOl by noting that as Ak in- 
creases {k2 increases) the corresponding transmission co- 
efficient T22 grows reducing the state uncertainty by the 
argument given for the resonances. The introduction of 
mixing terms, involved in 712,721, change the scenarios 
described above only quantitatively. The resonances will 
shift positions and the envelope of the concurrence as a 
function of Ak will now be non-monotone. As 77 depends 
on the determinant and the trace operations, one can di- 
agonalize the new 7 matrix and formally use equivalent 
expressions to the ones above. 

Although it is not the intent of the paper to propose 
a practical experimental setup to produce entanglement, 
resonance effects are ubiquitous for any quantum dot sys- 
tem coupled to external leads. The width of the res- 
onances can be controlled by the coupling of the dot to 
the leads and the resonance position in energy can be ad- 
justed, relative to the Fermi levels in the leads, by a gate 
voltage. The results of our paper show that the manipu- 
lation of the resonances will lead to the control of the de- 



gree of the entanglement. The realization of two separate 
electron channels in the same region has been addressed 
differently in the literature by tapping ^nto edge states 
that can provide two quantum numbers [8] for the incom- 
ing electrons. A gated quantum dot can be placed in 
the vicinity of the edge states with a controlled coupling 
so as to modulate the resonance characteristics of the 
dot. The outgoing electrons can be detected by gate elec- 
trodes that mix channels appropriately so as to change 
the measuring eigenbasis and compute for example the 
Bell inequalities or other measures of quantum correla- 
tions. This setup for the case of a spectrally structureless 
beam splitter has been described previously on the basis 
of current-current correlations in ref.[6]. 

This work analized the role of resonances in entangle- 
ment production in a quasi one dimensional two elec- 
tron system inspired in the electron-hole entangler of 
ref.[8]. Although we have restricted to elastic scatter- 
ing (no channel mixing), we found that by tuning the 
channel momenta or, equivalently, the resonant levels of 
a double barrier (through for example gate voltages) the 
quantum correlations associated with the post-selected 
scattering process can be manipulated in a controlled 
fashion. Needless to say, the response to resonances of 
electron-electron correlations is relevant in the context of 
electron hole entanglement Q . 
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